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Abstract 

Considering two-body integral equations we show how they can be dimcnsionally reduced by integrating 
exactly over the azimuthal angle of the intermediate momentum. Numerical solution of the resulting equation 
is feasible without employing a partial-wave expansion. We illustrate this procedure for the Bethe-Salpeter 
equation for pion-nucleon scattering and give explicit details for the one-nucleon-exchange term in the 
potential. Finally, we show how this method can be applied to pion photoproduction from the nucleon with 
TtN rescattering being treated so as to maintain unitarity to first order in the electromagnetic coupling. The 
^ ' procedure for removing the azimuthal angle dependence becomes increasingly complex as the spin of the 

. particles involved increases. 
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I. INTRODUCTION 



In cases when solving the Lippmann-Schwinger or Bethe-Salpeter type of equation is numerically 
involved, one often resorts to a partial-wave decomposition (PWD) in the center-of-mass (CM) 
frame. In doing so one can exploit the spherical symmetry of the interaction and perform the 
integration over the two-dimensional solid angle of the intermediate momentum analytically. While 
this reduces the equation's dimension by two, one has to deal with summing the partial- wave series, 
and hence this procedure is beneficial when only a few partial waves dominate. In the case when 
many partial waves must be taken into account, when restriction to the CM frame is not desirable, 
or when the potential is not spherically-symmetric, the partial-wave expansion is not helpful and 
one has to face the complexity of three- or four-dimensional integral equations. 

Fortunately, as had been noted by Glockle and collaborators [1, 2] in the context of the nucleon- 
nucleon (NN) interaction, the dependence on the intermediate momentum azimuthal angle factor- 
izes and can still be performed analytically without employing any kind of expansion or truncation. 
While this procedure has been successfully applied a number of times to the NN situation [2-4], 
here we would like to examine general conditions which potentials must satisfy to factorize the az- 
imuthal integration. We then apply it to solve a specific example of relativistic potential scattering 
in the pion-nucleon (irN) system and compare with the usual method of using the partial-wave 
expansion. 

In Section II we give the general requirements on the potential that allow one to remove 
the azimuthal angle dependence in the integral equation. In Section III we focus on the Bethe- 
Salpeter equation for 7rN scattering with one-nucleon-exchange potential and show in detail how 
the azimuthal-angle dependence can be integrated out in this case. Furthermore, in Section IV, 
we solve the resulting equation using a quasipotential approximation and compare the solution to 
the one obtained using the partial-wave expansion. In Section V we examine an extension of this 
approach to the calculation of pion electroproduction from the nucleon including the ttN final state 
interaction. Our conclusions are summarized in Section VI. 

II. CONDITIONS FOR EXACT INTEGRATION OVER THE AZIMUTHAL ANGLE 

The starting point in calculating observables of a two-body scattering process is an equation 
for the scattering amplitude. We shall assume relativistic scattering, in which case the equation is 
a 4-dimensional integral equation of the Bethe-Salpeter type: 



where T is the sought T-matrix, G is the two-particle propagator, and V is the two-particle- 
irreducible potential. Moreover, throughout the paper, q, q", q' stand for the relative 4-momenta 
of the incoming/intemediate/out going channel while P = p + k = p' + kl = p" + k" is the total 
4-momentum with k, k" , k! and p, p", p' the incoming/intermediate/outgoing momenta of particle 
one and particle two, respectively. 

In order to investigate the conditions under which the above equation can be integrated over 
the intermediate azimuthal angle we work in the helicity basis and only display the dependence on 
the azimuthal angle and helicity: 



T(q',q;P) = V(q',q;P)+i 




P)G(q";P)T(q",q;P) 



(1) 



T x , x (<f?,<p) = Vx.x(<f/,<p) + Y l I 




^V x , x »(<p',<p'')G(<ff')T x ,, x (<f/',<p). 



(2) 
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An important point here is that the two-particle propagator G can always be made independent 
of the intermediate angle ip" by choosing the total three-momentum along the z-axis, i.e., choosing 
the co-linear frame: P = (Po,Q, 0, P3). Furthermore, we shall observe that in the case when only 
spin-0 and spin-1/2 particles are involved, the azimuthal-angle dependence of the fully off-shell 
potential [13] in the co-linear frame is given as follows 

V yx ( V ', <p) = e"* A V v x . x {y' - V ) e^ (3) 

where A and A' stand for the combined helicities of the initial and final state, respectively. The 
half-off-shell potential then takes a very simple form: 

I half off shell A AW half-off-shell W 

where A is the helicity of the on-shell state. 

It is in this case, when conditions (3) and (4) are met, the exact integration over the azimuthal- 
angle can readily be done. First, by using (3) in Eq. (2), we see that the azimuthal dependence of 
the t-matrix is given by: 

T\>\((p' , ip) = e- iX '*' t yx (ip' - if) e 1 ^. (5) 
Since v and t only depend on difference tp 1 — ip, we expand them in a simple Fourier series: 

m m 
It is straightforward to show that their Fourier transforms, 

2tt 2tt 



/ g ^) *S = *A'A e"^ , (7) 



^A^-S+E-S^tS. (8) 



,(m) _ 



satisfy the following equation which does not involve the (^-integration: 

A" 



In principle, m runs to infinity and so we have an infinite number of equations to solve even 
though they are not coupled. Fortunately, since only the half off-shell potential is needed to solve 
the equations and it obeys condition (4), the corresponding Fourier transform is non- vanishing 
only for m = —A: 



(m) 
V X>X 



half-off-shell =^TOW0)| h „„ shell . (9) 



The scalar system is the simplest one where this procedure can be demonstrated. In that case 
the potential is a scalar function of scalar products of relevant 4-momenta: 



V(q', q; P) = V(q -q',P-q,P- q', q\q", P 2 ) (10) 

Given q = (qo, |q| sin# cos <p, |q| sin# simp, |q| cos#) and similarly for q' , we easily convince our- 
selves that, in the co-linear frame, the azimuthal dependence enters only through the product: 

Q " Q 1 = Q0Q0 — |q| |q'| [cos 6 cos 9' + sin sin 9' cos {ip' — ip)] (11) 
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FIG. 1: Diagrammatic form of a relativistic two-body scattering equation. 



and hence it is of the necessary form given in Eq. (3). Furthermore, in the half-off-shell case the 
momentum of the on-shell state, say q, can always be chosen along the z-axis, i.e., such that 9 = 0. 
Hence the half-off-shell potential is independent of azimuthal angles which fulfills condition (4) for 
the spinless case. The two-particle propagator G(q; P) = G(P ■ q, q 2 , P 2 ) is of course independent 
of (p in the co-linear frame. 

Once we have found that conditions (3) and (4) are satisfied, while G is independent of ip, the 
integration over <p can be done immediately. We will now show this more explicitly for the more 
complicated case of a scalar-spinor system. 



III. SPIN COMPLICATIONS: THE irN SYSTEM 

Consider the Bethe-Salpeter equation for the case of elastic scattering of a scalar with mass 
m w — the "pion" — on a spinor with mass tun — the "nucleon". We attribute the momenta 
p, p' to the nucleon and k, k' to the pion. The relative 4-momentum of the incoming channel is 
conveniently defined by q = ftp — ak, where Lorentz scalars a and ft are given by 

a = p ■ P/s = (s + m 2 N — m n )/2s , 

ft = k-P/s = (s-m 2 N + m w )/2s , (12) 

with s = P 2 . Similarly one defines q' = ftp' — ak' and q" = ftp" — ak" as the relative 4- momenta 
of the outgoing and intermediate state, respectively. In terms of these variables, the two-body 7rN 
Green's function of Eq. (1) is: 

r(n . P] = 1 (aP + g) • 7 + m N 

l? ' J (ftP -q) 2 - ml + ie(aP + qy-m 2 N + ie' [ } 

Projecting the equation onto the basis of the nucleon helicity spinors (defined in Appendix A), we 
obtain 

T^(q', q; P) = V&tf, q; P) + i £ / ^vg$W, q"; P) G^"\q"; P) T^q" , q; P), (14) 

A"p" 

where the helicity amplitudes are defined as 

T(M, q, P) = (1/47T) 4'\aP + q') T(q', q, P) u[ p \aP + q), (15) 
and analogously for V, while the defining equation for G^ is 

(aP + q) 7° G(q- P) 7 ° u[ p) («P + q) = 5 yx 5 p , p G^ (q; P), (16) 
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and hence 



Q(±)(q-p} = 

q + a^/s ± {E aP+q - ie) (/3y/s - q ) 2 - ^ p _ g + ie ' 



with E q = ^/q 2 + and u q = ^/ q 2 + m 2 . 

The most general Lorentz structure of the fully off-shell potential in the helicity basis can be 
written in the form [14] 

V^ // (q',q;P) = u(,(aP + c l ') [a(> + Ag" 7 ° + {A^ + ^V) 7 ■ p] u$(aP + q), (18) 

where Aj are scalar functions of the dot-products of the relevant momenta, i.e., 

A, = A i (q-q',P-q,P-q > ,q 2 ,q ,2 ,P 2 ). (19) 

Considering the dependence of these functions on the azimuthal angles of q and q', we see that 
- in the co-linear frame — it is given by the difference <// — (p, for the reason described below 
Eq. (10). 

The rest of the 93-dependence resides in the nucleon spinors. According to Eq. (18), in the 
co-linear frame we need to consider only Xv (©', ip') Xa(©> V 9 ) an d xt'(®'> V 5 ') °"3 Xa(©j V 9 ) where 
x's are the Pauli spinors (cf. Appendix A), 0, <p and 0', ip' define the orientation of ctP + q and 
aP + q', respectively. Since, 

x A ,(©', <p') XA (e, V ) = e - iA V' d v 2 „ (@ ') 4£(e) e ^"^'-^] (20) 

we observe that the (/^-dependence of these elements is of the desired form Eq. (3). And for the 
half-off-shell situation, where we can choose 9 = (hence = 0, in the co-linear frame) and use 
^AA"(0) = ^AA"> w e find the form, 

xt(eV)xA(0,¥>) = e-W-Wd)!,l(&), (22) 
xU 0, ^')^3Xa(O,^ = e-* x '-W (-l)V2-A d V;(e'), (23) 

which obeys the necessary half-shell condition Eq. (4). 

Therefore, we have demonstrated that the azimuthal-angle dependence of a pion-nucleon poten- 
tial in the co-linear frame always satisfies conditions (3) and (4). It is also apparent from Eq. (17) 
that the two-particle Green's function does not have any azimuthal dependence in that frame. 
Thus the integration over tp can exactly be done in the Bethe-Salpeter equation for ttN system by 
means of the procedure of Sec. II. 

Similar arguments apply in the case when both particles have spin 1/2, e.g., the nucleon-nucleon 
(NN) scattering. It should only be noted that in this case the potential satisfies conditions (3) and 
(4) with A = Ai — A2, A' = X[ — \' 2 . In other words, helicities of the two particles must be combined. 



IV. NUMERICAL RESULTS 



The standard route to solution of a potential scattering equation such as Eq. (14) is to decompose 
it into an infinite set of equations for partial- wave amplitudes, see e.g. [8, 9]. The advantage of doing 
a partial wave decomposition is that the equation for each partial wave is of 2 lesser dimensions 



5 



k p' 



k' 



FIG. 2: One-nucleon-exchange irN potential. 



than the original equation, while the partial-wave series is usually rapidly converging, hence only 
the first few partial-wave amplitudes need to be solved for. 

On the other hand, solving for the full amplitude directly has its own important benefits. And 
if the exact azimuthal-angle integration can be done a priori, the numerical feasibility of this 
approach becomes comparable to the PWD method. 

In this section we would like to compare the two methods for the example of solving a relativistic 
equation for the irN system. For our toy-calculation potential we take the one-nucleon exchange, 
Fig. 2, and use the instantaneous approximation, thus neglecting retardation effects in the potential. 
The latter approximation allows us to perform the relative-energy (go) integration such that we 
are left with a relativistic 3-dimensional Salpeter equation: 



r A ^(q',q;P) = ^(q',q;P) + W ^^(q', q"; P) Gg? (q"; P) rft(q", Q! P), (24) 

\" nil J 



\"p" 

where the equal-time two-particle propagator in the CM system is given by 



4^(|q|; v^) = 2, r ^GW( ffi P) = — rTlh M < 25 ) 

J-oo 2vr u q {-py/s + E q +u) q - ic) 

This 3-dimensional equation for ttN has been described in detail and solved using a PWD in 
the CM system by Pascalutsa and Tjon [5-7]. We, on the other hand, solve this equation by 
using the framework of the two previous sections to reduce the 93-integration analytically and solve 
numerically the resulting 2-dimensional integral equation for the m-th Fourier component of the 
full amplitude: 

t^' p ^\,6\H,6)=v^ p ' p (\ (l '\,e',\c l \,e) 

d\q" 



q "| 2 f de"v^f p '(\a'\, e', \d"\,9") G^\\c l "\)t^{ p '' p (\c i %e", \q\,e), (26) 

J 



where, without loss of generality, we have also assumed the CM frame. The explicit form of the 
Fourier transform of the one-nucleon-exchange potential is worked out in Appendix B. 

Let us emphasize that it is necessary to solve for only one of the Fourier components (either 
m = —1/2 or m = 1/2), the other ones either vanish or can be obtained by relations due to the 
parity and time-reversal invariance. 

We solve Eq. (26) by the Pade approximants as in Refs. [6, 7] thus maintaining exact elastic uni- 
tarity. The numerical integrations are performed by the Gauss-Legendre method. The integral over 
|q"| in Eq. (26) contains the cut singularity at |q"| = ^[s — (mjy — m,,-) 2 ]^ — (mjv + m n ) 2 ]/4:S = q, 
which is handled by the well-known identity: 



L 



iH fM^ =v r M \l^i- mm (27) 
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where V denoted the principal-value integral. When computing the latter the integration region is 
divided into two intervals: |q| G [0, 2q], and |q| £ (2q, oo). The Gaussian points are then distributed 
separately for each interval to make use of the property that an even number of Gaussian points 
falls symmetrically with respect to the middle of the interval hence the singularity in the middle of 
the first interval is avoided. The polar angle integration is straightforward for both the principal 
value term and the imaginary contribution. We find it sufficient to use 16 Gauss points for the 
momentum integration and 8 points for the polar-angle integration. Upon increasing the number 
of points to 32 and 16 respectively, the results change by less than 0.5% in the considered energy 
range. In all cases we found that 6 iterations combined with the use of Pade approximants works 
extremely well. 

After we solve Eq. (26) to find the full ttN T-matrix, we can of course also find the partial wave 
amplitudes: 

T£ V (|q'|,|q|) = r^r A ^(|q'|,|q|^)4 A (0), (28) 
Jo 

where 6 is the angle between q and q'. We then investigate the convergence of the partial wave 
series: 

T^(|q'|,|q|^)=^(j+I) T&» (|q'|, |q|) 4aW • (29) 
J 

In particular, in Fig. 3 and Fig. 4 we plot the on-shell values of |T"w?| 2 compared with the truncation 
of the partial- wave series for 3 terms and 5 terms (i.e., J = ^, . . . , | and J = \, . . . , § respectively). 




FIG. 3: Angular dependence for \Ttt\ 2 at E^ AB = 300 MeV. Solid line is the full calculation, dashed and 

2 2 

dotted are the resumming of partial terms. 

In order to compare the computational efficiency of the two methods, we compare the number 
of partial waves needed to achieve convergence in the PWD method with the number of Gauss 
points for the polar-angle integration which appear in the "w/o PWD" method. 

The figures show that the effect of truncations of the partial-wave series increases with the 
angle (Fig. 3) the energy of the incoming tt (Fig. 4). In our particular case of one-nucleon exchange 
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[MeV] 



FIG. 4: Energy dependence for \T+ti\ 2 at 9% M = n. The lines arc defined the same as in Fig. 3 
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computing 5 or more partial wave amplitudes is sufficient to reproduce the full result to a 1 per 
cent accuracy in a broad energy domain. Thus, in this case, the efficiency of the two methods is 
comparable since we need 5 multipoles versus 8 Gauss points of the polar-angle integration. 

It is important to emphasize that the ability to do the azimuthal-angle integration analytically 
is necessary to achieve comparable efficiency. We have checked that it usually takes at least 16 
Gaussian points for the azimuthal integration which slows down the calculation by more than an 
order of magnitude. 



V. EXTENSION TO PION PHOTOPRODUCTION 

Our procedure for performing the analytic 92-integration is applicable in the photo- or electro- 
meson production to first order in the electromagnetic coupling. Here we describe the extension 
to the case of it photoproduction within a simple final-state-interaction model [6, 10]. The model 
begins with the following coupled channel equation: 

(T n7T ^~7T7 j J ^7T7T ^7T7 \ , / ^7T7T ^7T7 j ( G \ / T n7T T-^y 

Ty n J \ V~fir ^77 / \ ^7"7r ^77 / \ G-y J \ T^-k -^77 

where T and V are the amplitudes and driving potentials of the 7rN scattering (inr), pion photo- 
production (771"), absorption (ttj), and the nucleon Compton effect (77), respectively. The above 
equations are solved up to first order in the electromagnetic coupling e, hence preserving two body 
unitarity to this order only. 

In solving the photoproduction scattering equation we calculate first V n7T as described for 7rN 
scattering and we then iterate in the following manner: 

where we used = T 77r from time-reversal invariance. 

This solution procedure is obviously suitable for our case since the half-shell V n7 has a simple 
azimuthal angle dependence similar to the case of (see Eq. (B7)). The reduced kernel (see 
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Eq. (Bll)) has two terms rather than the one term in the 7rN case due to the "complication" of 
having to couple a spin 1 photon to spin 1/2 as opposed to coupling a spin meson to spin 1/2. 
For example, if one considers the nucleon u-channel exchange (compare to the 7rN case in Eq. (B2)) 
the half-shell photoproduction potential can be written as: 

V$LW, Q) = iVfa{<&, |q', Qo, H,9')e-^'-^' + 2 V^ a (q' , |q', q , |q|, 9') e -*(*+W (32) 

where a = ±1 represents the helicity of the incoming photon. 

One sees that when Eq. (32) is iterated in Eq. (31) two de-coupled scattering equations are 
obtained (each corresponding to \V or 2V). For each of these equations, one can show that the 
corresponding <p' dependence re-appears after doing the tp" integration and therefore once again 
we can perform the azimuthal-angle integration analytically. As in the 7rN case, the resulting 
"reduced" kernels obey 2-D integral equations. 

As a check of our procedures we calculated the u-channel contribution to pion photoproduction 
using the analytic azimuthal-angle integration along with 2-D numerical integration and compared 
to the results of Refs. [6, 10] obtained using the multipole expansion. At of 300 MeV with five 
multipoles we found agreement to better than 1% over a wide angular range. 



VI. CONCLUSION 

In recent years Glockle and collaborators [1, 2] introduced a method which greatly simplifies 
the numerical integration of two-body scattering equations without performing the partial-wave 
expansion. The method exploits a certain azimuthal symmetry of the potential thus allowing exact 
integration of the azimuthal dependence. In this paper we have established general form of the 
azimuthal-dependence of the kernel which allows for this procedure to go through. We have argued 
that these conditions is in general applicable to any system of spin- and/or spin- ^ particles 

We have applied this method to the case of 7rN system and . With some extra effort it can 
be applied to higher spin systems, however the procedure becomed increasingly complex with the 
increase of the spin of the involved particles. We have successfully applied the method to pion 
photo- and electro-production from the nucleon, however only to the leading order in electromag- 
netic coupling. 

Even though we have used the Salpeter equation for or numerical exercises, the method can 
of course be applied to the full 4-D Bethe-Salpeter equation, which for the ttN system has so far 
been solved in partial waves only [11, 12]. Performing the azimuthal-angle integration analytically 
greatly facilitates finding the full solution and makes the numerical feasibility of this approach 
comparable to tinging the solution using the partial-wave expansion. 



APPENDIX A: HELICITY SPINORS 

We define the four-component nucleon helicity spinors as follows: 

^• p) =7kU^P^) ^> (A1) 

where A = ±1/2 is the helicity, E p = \J\v\ 2 + m2 N is the energy, 6 and ip are the spherical angles of 
the 3-momentum p, and xa is the two-component Pauli spinor. The positive- and negative-energy 
nucleon spinors in the convention of Kubis [9] are defined as follows: 

u ( ±\p)=u x (±E p ,p), (A2) 
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They satisfy the following orthogonality and completeness conditions: 

4 W (p)«^ ) (p) = 6 pp ,5 xy , (A3) 
X;«? ) (P)«I W (P) = 1- (A4) 

The Pauli spinors along the z-axis are given by 

Xl/2(0)=(J), X-1/2(0)=(J), 

while along an arbitrary direction 9, ip they can be obtained using the Wigner rotation functions: 

X a(^)=e/;?w^ m> »(o), 

or, explicitly, 

/ cos(#/2) \ /-e-^sin(0/2)\ 
X^OM = U^in(0/2) ) ' X-i/2(f,v) = [ cos(e/2) ) ■ 

APPENDIX B: AZIMUTHAL DEPENDENCE OF ONE NUCLEON EXCHANGE 

As an example, we consider the u-channel nucleon exchange potential given by the graph in 
Fig. 2 and the following expression, 

v«, ,; p) = fa 7 ■ vp - ^v ( ,r^ P c + r! l r"'i + : w -' 5 -' ■ ««• - «»• < m > 

4m^ [(a - (3)P + q + q'\ 2 - m z N + is 

where a and f3 are defined in Eq. (12). For simplicity we choose the CM frame, where the potential 
in the helicity basis takes the form: 

~ 2 



■l>'l> i,J '■' , 9-nNN 1 

lGirmjj u — m 2 N 



xu}tf)[Mr 1 + M$> 7°]«^(q ff ) (B2) 



with 



M? p "(p',p") = m N [2j7s{p' Q +p'>) - s - 2p'.p" 

+p >2 + p "2 + m 2 N + (p / _ ^^(pg _ p » Ep „) 

-v^(Po " + PO " / (B3) 

M^'(p',p") = v^[2v^(Po +Po) - s - 2p'-p" 
-p' 2 - P" 2 ~ 3m 2 N - (p' - p'E p ,)(pZ - p"E p „)} 

+ {p" 2 +rn 2 M-p'E p ,). (B4) 

The azimuthal dependence cj>" arises from Dirac spinors and from various scalar products in- 
volving the the four vector q". Choosing the vector part of the total momentum P to be along the 
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z-axis (or to be zero in the CM frame) allows the </>" dependence, for the fully off-shell potential, 
to be displayed in the form: 



■2 



^ ^iK, |q'U', |q"|,^, O + d 2 (q' , |q'|, |q"|, 0', 0") cos(0" - 0') ' ( j 

where 7J^„ n , N p > p >i, and are factors which depend on the type of the diagram and of the 
exchanged particle, but are independent of the azimuthal angle. The quantities, 

m=-l/2 



Ng'g" = y (-Eg' + rn N ){E ql , + m N )/AE ql E qll 

are factors which result from the helicity spinors. In Eq. (B5) we have employed the usual trigono- 
metric relation between two arbitrary directions defined by q" and q': 

cos e q » q / = cos 6" cos 9' + sin 6" sin 9' cos(<£" - </>'). (B6) 

It is easy to see that the fully off-shell potential in Eq. (B5) has the azimuthal dependence of 
Eq. (3). Furthermore, in iterating Eq. (14) the quantization axis is defined by the on-shell relative 
momentum q (i.e. 9 = 0), hence Eq. (B6) reduces to cos6 q / iq = cos 9', therefore the half- off- shell 
potential reduces to: 

Vfctf, q) = v(,{(q' , |q'Uo, |q|, ^Je"^'"^ (B7) 

which if of the form of the result in Eq. (4). Therefore, the azimuthal angle dependence can be 
removed from the Bethe-Salpeter equation for this case. We achieved this result by explicitly 
displaying the azimuthal angle dependence and align P with the z-axis so that only 7 3 and 7 
appear in V£,£. The presence of 7 1 or would introduce additional azimuthal angle dependence 
in the spinor matrix elements and make the algebra much more complicated. 
For the -u-channel nucleon exchange the coefficients d{ are: 

di (p',p") = p' 2 + p" 2 + s - 2V~s(p' + p'o') + 2p' p / ' 

-2|p'||p"| cos 0' cos 6>"-m^ (B8) 



d 2 (p',p") = -2|p'||p"| sin 9' sm9" (B9) 

From these relations one can exactly identify the angular dependence of potential given in Eq. (Bl) 
in the 4-product, 

/ // / // in 
P-P =PoPo -P P 

= PoPo- |p , ||p"|cos0 , cos0" 

-|p'| |p"| sin 9' sin 9" cos(e/>' - <f>"). (BIO) 

The relative momenta q' and q", defined in SECTION II, are to be introduced in Eq. (B3) - 
Eq. (BIO) by p' = aP + q' and p" = aP + q" 
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where 

m=-l/2n=0 

After applying standard trigonometric manipulations: 



1 I 2 N r2n j n 

x, M \M'\y,#)= E E«(^')j[ d + 7Js/ ix ' m)x - (B11) 



cos 2 ™- 1 



and 



cos 2n 



22n-2 



1 / 2n 



cos(2/;- 1)0+ ( 1 J cos(2n-3)0 + ...+ ( 2 ^_ ^ leosfl 



+ 



1 



2 2n \ n J 2 2n - 1 



cos2/i0+( 2 ™ J cos(2n-2)0 + ... + ( 2 _^ ) cos 20 



the integral over the azimuthal angle of the intermediate momentum in Eq. (Bll) can be reduced 
to integrals of the following type: 



f 2n #( 
' n ~Jo ~ 



d(j) cos(m(f))e 1 ' 



d(j) cos(m</>) cos(n(; 



7o i 



+ a cos Jo 1 + a cos < 

For values \a\ < 1, this definite integral can be evaluated analytically to obtain: 



l m,n 



IT 
J) 



6-1 



m+n 



+ 



6-1 



\m— n\ 



(B12) 



(B13) 



where 6 = \/I — a 2 . 

The results given above in Eq. (Bll) work for all standard particle exchanges in the s, t, or u 
channels. Furthermore, it should be noted that additional azimuthal angle dependences introduced 
by various form factors can easily be handled by simple algebraic methods. The maximum power of 
cosx needed for a particular diagram may increase (for example, N=2 for u-channel A exchange). 
In addition, v will, in general, contain a sum of various terms corresponding to each diagram 
included. However, all of these terms can be evaluated using Eq. (B13). In addition as noted 
earlier, this procedure is not at all affected by the equal time approximation and can be applied in 
the same manner to the full 4-D Bethe Salpeter equation. 
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